Quantum relaxation in open chaotic systems 
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Using the super symmetry techniqu e, we analytically derive the recent result of Casati, Maspero 
and Shepelyansky [ cond-mat /970610q ] according to which the quantum dynamics of open chaotic 
systems follows the classical decay up to a new quantum relaxation time scale t q ~ y/t c ts- This 
scale is larger than the classical escape time t c but still much smaller than the Heisenberg time tn- 
For systems with orthogonal or unitary symmetry the quantum decay is slower than the classical 
one while for the symplectic case there is an intermediate regime in which the quantum decay is 
slightly faster. 



PACS numbers: 05.45. +b, 05.60,+w, 72.20.Dp 



The classical decay probability of a generic weakly 
open chaotic system obeys the exponential distribution 
P c i(t) cx e~*/* c where the mean escape time t c character- 
izes the effective coupling to the outside. Motivated by 
recent experiments on mesoscopic cavities or microwave 
billiards there has been renewed interest in the problem 
of quantum life times For example the quantum 

properties of "chaotic" maps with absorption were 
investigated. Recently, also analytical results for the sta- 
tistical distributions of the complex poles of the scatter- 
ing matrix |^,|| or of the eigenvalues of the Wigner-Smith 
matrix of time delays || were found for the generic prob- 
lem of chaotic scattering. 

The problem of current relaxation in disordered met- 
als which is similar to the quantum decay of an initially 
localized wave packet inside a chaotic cavity or a disor- 
dered region has been investigated by different analytical 
approaches |9|-[ll| in the framework of the nonlinear a 
model (replica or supersymmetry variant Jl2[). In these 
works it was shown that the classical exponential decay 
is strongly supressed by quantum effects for time scales 
larger than the Heisenberg time tn giving rise to a log- 
normal distribution of relaxation times. Muzykantskii 
and Khmelnitskii demonstrated that this is due to 
a nontrivial saddle point of the er-model. For the case 
of open one-dimensional geometries they obtained for 
t > tg the behavior P(t) ~ exp[— g ln 2 (i/t#)] where 
g = tH/t c ^> 1 is the conductance (in units of e 2 /h). 
This result has an important relation to the probability 
to find a 'nearly localized state' in a normally metallic 
sample |fTof . Also the quantum time evolution of open 
chaotic cavities was studied ]T^ | giving a power law de- 
cay for t^$> tjj/ min {Tj} where < Tj < 1 are the trans- 
mission coefficients of the barrier by which the cavity is 
coupled to the outside. 

Recently, Casati, Maspero and Shepelyansky 0] sur- 
prisingly found that for a quantum kicked rotator model 
with absorption Q significant deviations from the clas- 
sical behavior already appear at an earlier quantum re- 
laxation time scale t q ~ y/t c t^ -C tf[. Their argument 
jjj is based on the complex eigenvalues of the non uni- 
tary time evolution operator being typically distributed 



in a narrow ring of width E c = 1 jt c M inside the unit 
circle. (t c 3> 1 is measured in units of the kick period.) 
Then t q can be identified |Q as the inverse of their typi- 
cal distance in the complex plane. This picture is indeed 
supported by numerical quantum simulations [0 for the 
kicked rotator. 

In this work, we present analytical results for a similar 
model by mapping it onto the supersymmetric nonlinear 
cr-model |l^] which is possible due to a recent progress of 
Altland and Zirnbauer |L4| for this type of systems. The 
cr-model in the zero-dimensional limit also applies to the 
case of a chaotic cavity coupled to external leads [jl5| . We 
clearly confirm the findings of Casati et al. that the new 
time scale t q is indeed highly relevant for the problem 
and, additionally, we find at t ~ t q qualitatively different 
quantum effects for the three symmetry classes of random 
matrix theory |l6| which are characterized by the index 
(3 = 1 for the orthogonal case (systems with time rever- 
sal symmetry and no spin mixing), (3 = 2 for the unitary 
case (broken time reversal symmetry) and (3 = 4 for the 
symplectic case (time reversal symmetry and strong spin 
mixing). This result supports the interpretation that the 
effect of weak localization (or anti-localization for (3 = 4) 
can also be observed in open systems with absorption. 

We consider the quantum dynamics \ip(t + 1) > = 
S\ip(t) > of a generalized random phase kicked rotator 
model with the time evolution operator introduced in 
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where I, I are the quantum numbers of "angular momen- 
tum" beeing conjugated to the angle 9. As in Refs. [^]|l]] 
the Z-space is finite: —N/2 < 1,1 < N/2 introducing 
effective absorption at the boundaries. We consider ran- 
dom phases /i; jlTj and a quite general periodic kick po- 
tential V(9) (with a finite number of harmonics). The 
different symmetry classes are encoded in the symme- 
tries of V{6) @, i.e. V{&) = V{-9) for (3 = 1 and 
V{9) = V (9)t 2 + J2l=i°vV u (9) for = 4. Here a v 
are the Pauli matrices and V v {6) is even (or odd) for 
v = (or v = 1,2,3). In the following, we consider 
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the phase averaged quantity P(t) = (|< 0| S* |0 >| ) M to 
describe the decay of a quantum state initially localized 
at the site 1 = 0. For short time s cale s this probabil- 
ity decays diffusively as P(i) oc 1/y/Dt (with diffusion 
constant D = (V'(9) 2 )g ^> f) whereas for longer time 
scales and large system size (t, N ^> D) quantum lo- 
calization leads to the saturation P(t) oc l/£ with the 
localization length £ = (3D/2. Here we concentrate on 
the case of a system size N being much smaller than £ 
(i.e. t c ~ N 2 /D <^C A" = f#) and on time scales t > t c . 

We first present and discuss our main results before we 
outline some basic steps of the approach. We find that 
the first quantum corrections for t c < t < t 
be cast in the form 
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where for the Od limit (k as N/2 or for a chaotic cav- 
ity) we have introduced the "generalized conductance" 
moments by g u = iy, ^ = f , 2, . . .. Here the trans- 
mission eigenvalues < Tj < 1 describe the effective 
coupling strength of the cavity with the boundary. For 
the fd limit, g\ = g 2 = tt 2 D/2N is (up to a numerical 
factor) the classical conductance from the site to the 
boundary. The universal functions Cp(u) have the form 
(inset of Fig. [TJ) 
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Cl (u) = l d X j Q dy {1 *f +x2y2)2 « 

x exp [u (2x - (1 - x 2 + x 2 y 2 ))] 
= l + u 
C2(u) = sinh(w)/w 
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C 4 (u) = Ci(-u/2) = l-in+lr« 2 
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For /3 = 1,2 the quantum probability P(i) is above its 
classical value Pei(i)- The criterion ln[P(t q ) / P c \(t q )} = 
0.1 (sec Q) to define the quantum relaxation time scale 
t q leads to t q w 0A5y/N/E c ^ for /? = 1 and f 9 w 
1.24y/N/E C:2 for /3 = 2. The numerical factor for /3 = 1 
is indeed close to 0.38 found in [jlj. We note that for 
/3 = 2 the function C2(w) has only a quadratic correc- 
tion for small u. The situation for (3 — 4 is particularly 
intriguing because here the quantum probability is ini- 
tially even below the classical value. The function C&(u) 
has at u m in ~ 3.03 (tq, m i n ~ 2.46-y/ 'N / '-Be, 2) its minimum 
value 0.488 and it crosses the classical value 1 again at 
w cr w 7.36 (tg, cr ~ 3.84-\/A r /£ , c ,2). ft seems that the lin- 
ear term in the function Cp(u) can be viewed as a weak- 
localization correction (anti-localization for = 4). In 
Fig. [j], we also show for the Od case with Tj = 1 the 
accurate result which is given by more complicated inte- 
grals (see below). 
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1. Logarithm of P(t) for the three symmetry classes 



0. 



Tj = 1. The full lines are obtained from 



1, N = 2000 and all transmission eigenvalues 
0) for (3 = 2 or the 
corresponding integrals for = 1,4. The dashed line shows 
the classical exponential decay and the two dotted lines for 
P = 1,2 correspond to Eq. (|2). For (3 = 4, Eq. (Q) coincides 
with the full line. The inset shows the functions Cp(u) given 
in Eqs. (|)-(|). 

To derive these results, we have applied the supersym- 
metric technique [|l2],[l5) which has recently been gener- 
alized 0,^9) to treat random phases instead of gaussian 
disorder. Repeating the steps described in Ref. |l9), we 
can express the Laplace transform P(u>) of P(t) as a func- 
tional integral of the type 



PQ/(Q(0)) c- £[c?1 . 



(6) 



Here the integration is done over a field of 8 x 8 super- 
matrices Q(l), —N/2 < I < N/2 with the non linear 
constraint Q 2 — 1 and particular symmetries for each 
universality class |Lq ). /(Q(0)) is a preexponential fac- 
tor that depends only on the Q-field at site 0. The action 
in (^|) has the form 



C[Q] = f Str 8JV \n(B(uj)+iQ) , 



B(oj) = iA 
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The number d = 1 (2) for j3 = 1, 2 ((3 = 4) measures 
the spin degeneracy and the supertrace extends over an 
8iV-dimensional super space. Q is an operator containing 
the Q(Z)-fields in its diagonal blocks and Uq is a matrix 
with elements < l\ e~~ iV< - 6 ^ \l > (8I4. The block structure 
in (||) refers to the grading for advanced and retarded 
Greens functions with the matrix A having the entries 
+1 (—1) in the upper (lower) diagonal block. As in [p^| , 
we expand the action in the limit of long wave lengths 
and long time scales which gives £[Q] ~ Cb[Q] + £id[Q] 
where 
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is the standard one-dimensional a model action. Here 
the supertrace without subscript acts on 8 x 8 superma- 
trices. The term £b[Q] which was absent in @Q arises 
from the boundary absorption because the operator XJq 
is not unitary due to the cutoff in I space. According to 
this we can write B(0) = B± + iAB 2 with hermitian ma- 
trices B\ and B 2 - Note that B2 docs not vanish because 
U$ is not unitary. The boundary part of the action is 
then determined by the eigenvalues < < 1 of the 

hermitian matrix = A" 1 / 2 AB 2 A~ 1 / 2 (with A = 
B\ + (1 + B 2 ) 2 ). These eigenvalues have the meaning of 
transparencies of coupling channels to the outside. Their 
precise distribution depends on microscopic details like 
system size and the particular choice of the kick poten- 
tial V(9). The eigenvectors with non vanishing T. have 
typically a support on the sites close to the boundary and 

the related boundary conductance — J2j Tj scales 
like the effective bandwidth of Uq: ~ \f~D. We have 
verified this behavior by a numerical evaluation of 
for the standard kicked rotator. Therefore we can write: 
C B [Q] ^L B (f(°\Q(^)) + L B (f^,Q(-f)) with 



Mf(°>,Q) = J £Str Ml 



iT^ AQ) 



(10) 



and AQ = j(QA + AQ) — 1. The sum runs over all 
non vanishing eigenvalues associated to one boundary. 
We note that for the S'-matrix approach of Refs. [ |15|j20| 

exactly the same action is obtained where are the 
transmission eigenvalues of a tunnel barrier which cou- 
ples a mesoscopic sample to an ideal quantum wire pof . 

The functional integral (|^) corresponds to a path inte- 
gral which can be evaluated by solving a diffusion equa- 
tion in Q-space pl|^2|j2^ ] . Therefore we rewrite (||) as 



P(oj)= / dQ f(Q) F 2 (Q,N/2) 



(11) 



where the function F(Q, I) 
differential equation |2l|,|lS 



is determined by the partial 



F (Q,l) 



J df p(f, I) expf-Li 



(T, Q)] 
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where p(T, I) is a probability distribution of transmis- 
sion eigenvalues Tj which fulfills a certain Fokker-Planck 
equation (known as DMPK-equation due to Dorokhov 
p4| . and Mello, Pereyra, Kumar [p5| ) with the initial 
condition p(f, 0) = 5(f-fW). p{f,l) describes the sta- 
tistical transport properties of a quasi one-dimensional 
disordered wire in series with a tunnel barrier with trans- 
parencies Tj°\ At first sight ( p"3| ) seems to be more 
complicated due to the increased number of integrations. 
However, in the metallic limit, we can expand (|l^) in 
powers of AQ with the self averaging transmission mo- 
ments gi,ff2,93> ■ ■ ■ as prefactors. Their "quantum" fluc- 
tuations are of order unity and have only an effect for 
t > tn- Therefore we can replace g v by their average 
values and omit the T-average. These g„-averages are in 
the classical limit determined by a set of differential equa- 
tions which can be derived from the DMPK-equation [^6| . 
To determine F(Q, I) for ui 7^ we use the expression for 
F (Q, I) as an ansatz where the g v are now parameters to 
be determined as a function of ui. The w-term only mod- 
ifies the equation for g\ giving: g[(l) — ~(2/D)g\ — iuu 
and g' 2 (l) = (4/-D)(<? 2 — 2g\g 2 ). Omiting the details, 
we mention that the explicit solutions determine F(Q, I) 
and thus provide a closed expression for P(lo) as one 
Q- integral (11). Using the standard parameterizations 
for Q introduced by Efetov |12|1 , we can express JIl| ) 
as an integral over two (/3 = 2) or three (f3 — 1,4) ra- 
dial parameters. We can perform the integrations for uj 
(from the Fourier transform) and for the effective vari- 
able s = Str(AQ) in a saddle point approximation which 
is justified for t 3> t c . Keeping the first two terms with 
gi and g 2 in F(Q, I) we obtain our main result (|^-||) for 
the Id case. The situation for the Od case is much easier, 
here we can simply insert the given 'boundary' transmis- 
sion eigenvalues and perform the w-integration. For lack 
of space, we only state the result for (3 — 2 
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1 — x To 
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and the initial condition F(Q,0) = exp[-L B (T(°>, Q)]. 
Here Aq denotes the Laplace operator in Q-space (with 
the precise notations of Ref. |2^]). The general solution 
of (|l^) for arbitrary frequencies is an involved mathemat- 
ical problem. First, we consider the solution Fq(Q,1) for 
the case uj = 0. For this, we note that exp[— Lb(T, Q)] as 
a function of Tj and Q exactly coincides with the gener- 
ating function (2.3) of Ref. |23|] which was used to prove 
the equivalence of the a model [ pl| , p2[p2] | and Fokker- 
Planck approach [p4|-|2lf| for quasi one-dimensional disor- 
dered wires. According to the argumentation presented 
in UJ, Fq(Q, I) is exactly given by 



The corresponding expressions for = 1, 4 have a sim- 
ilar structure with two integrations. The curves shown 
in Fig. ^ were obtained from a numerical evaluation of 
these integrals. They also lead to our principal result 
@"@ ^ we expand the logarithm in (JlJ) up to second 
order in T. The expansion parameter here is in principle 
t/N ~ t/tii <C 1. However, one can estimate that the 
third order term gives a contribution oc i 3 /{t c t 2 H ) which 
has to be smaller than unity because of the exponential 
in (0). Of course the same criterion holds for the Id 
case if we restrict ourselves to the first two moments g± 
and g 2 . 
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In summary, we have found that for open chaotic sys- 
tems the first quantum corrections to the classical re- 
laxation process appear at a quantum relaxation time 
scale t q ~ y/t c tn with different effects for each univer- 
sality class (Fig. |l|). This scale is determined by the 
second moment of transmissions eigenvalues Tj describ- 
ing the effective coupling strength of the initial site with 
the boundary. It would be very interesting to relate this 
finding more clearly to the physical mechanism suggested 
in Ref. jjj according to which t q is the time scale at 
which the quantum discreteness of the complex eigen- 
values exp(iEj — Tj/2) of the non unitary time evolution 
operator S || can be resolved. We emphasize that in 
view of the universal a model formulation our results ap- 
ply not only to the kicked rotator model (0) but also to 
chaotic cavities (corresponding to the zero-dimensional 
random matrix limit) and to quasi one-dimensional dis- 
ordered wires. In this case one should consider the time 
evolution of a wave packet of plane waves in an energy 
interval of size U/t where r is the elastic scattering time. 
The typical extension of a such a wave packet is just the 
mean free path which is in any case the smallest length 
scale that can be resolved by the standard a model [|2| . 

Due to the almost identical a model action it is im- 
portant to understand the relation of our results with 
the approach of Ref. |l(J where mainly the limit t > tjj 
was considered. A recent careful analysis 1271 of the sad- 
dle point approach pioneered in Ref. |fL0| indeed gives 
for the regime t q <C t <C ta the behavior \nP(t) w 
~{t/t c ) [1 - t/(j3gt c )] confirming Eqs. (|)-(|) for u > 1. 
Furthermore, for t > tu we can state that the log-normal 
behavior found in [p|-pr| should also apply to the average 
decay rate for the kicked rotator model. However, for 
very long time scales one should also focus on the distri- 
bution of the decay function because for a given sample 
the decay is then again exponential with a decay rate 
given by the minimal Tj Q . 

Concerning the zero dimensional limit, for (3 = 

1 is in principle also contained in the exact integral ex- 
pressions of fl3| ]. However, since the corresponding limit 
was not worked out there the time scale t q remained un- 
detected. We emphasize that here the Tj are given model 
parameters and E c ^ might parameterically be smaller 
than E c ,i if all Tj <C 1. We mention that very recently 
Savin and Sokolov p8| independently also found the time 
scale t q in the frame work of the supersymmetric ap- 
proach. Their results which apply for the Od case with 
unitary symmetry completely agree with our findings (|^) 
and @). 

The author acknowledges D. L. Shepelyansky and 
B. Georgeot for fruitful and inspiring discussions. 
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